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ABSTRACT 

An axial crack in a cylindrical shell is investigated by use of 
a 10th order shell theory, which accounts for transverse shear defor- 
mations as well as a special kind of orthotropy. The symmetric prob- 
lem is formulated in terms of two coupled singular integral equations, 
which are solved numerically. The asymptotic membrane and bending 
stress fields ahead of the crack are found to be self similar. Stress 
intensity factors are given as a function of the shell parameter for 
various values of the ratio crack length to shell thickness. Consider- 
able differences from 8th order shell theory results are found for the 
bending stresses, while the membrane stresses of the 8th order theory 
seems to be a lower limit reached for very thin shells. 


Part of the present work was carried out during a visit to Lehigh 
University supported by NSF under the Grant ENG 73-045053 AOl and 
NASA-Langley under the Grant NGR 309-007-011. 
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r. immcnoN 

Plate and shell theories based on Kirchhoff'.s assumption only en- 
ables the satisfaction of two plate-type boundary conditions. Broadly 
speaking the importance of this restriction depends on the variations 
of the field quantities compared to the plate or shell thickness. 

Thus serious shortcomings can be expected in problems where steep grad- 
ients are encountered, e.g., crack problems. A number of investiga- 
tions [1-4] concerned with a crack in an infinite plate in bending have 
revealed considerable differences between the 4th order bending theory 
solution and solutions obtained by use of 6th order Reissner-type 
bending theory [5]. The differences concern both the magnitudes of 
the stresses and their distribution around the crack. In the light 
of these results it seems to be of considerable interest to supplement 
existing results for cracks in shallow shells described by 8th order 
shell theory with calculations, which explicitly incorporate the effect 
of transverse shear. 

In the present paper a brief derivation of shallow shell field 
equations including transverse shear is given. The procedure follows 
closely that of Naghdi [6], but a certain type of orthotropy is in- 
cluded here. It is demonstrated, how these orthotropic equations can 
be obtained from the isotropic equations merely by use of suitable 
variable transformations. This also holds for the 8th order shell 
theory used in previous investigations [7-10], and these solutions 
can therefore be given a more general interpretation. 
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Although the shallow shell equations are fairly general, the rest 
of the paper is confined to an axial crack in a cylindrical shell. 

The method of solution is that of singular integral equations obtained 
by the complex Fourier transform. An investigation of the asymptotic 
stress field around the crack tips is included, and extensive numer- 
ical results are given. 


Z. FUWPAMEMTAL EQUATJOfJS 


In this section a brief derivation is given of the fundamental 
equations for a shallow elastic shell. A special kind of orthotropy 
is accounted for through the parameter 6. In the isotropic case the 
equations are those given by Naghdi [6]. The middle surface of the 
shell is described in the cartesian coordinate system of Fig. 1 as 
Z=Z(X^), a=l,2. The stress resultants are given in terms of their 
components in the {X„,Z} system. A different approach has been used 
by Sih and Hagendorf [11] for an isotropic spherical shell. 

When only a vertical load q{X ) is included, the equilibrium equa- 
tions are 




= 0 
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shear forces. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


are the membrane forces the moments and the transverse 


(2.4) 


For a shallow shell the strains are defined by 


„ 1 


'ctp 




where U and W arc the displacement components in the {X ,Z} system. 
The normals to the shell in its original configuration change direc- 
tions by the angles The slope of the middle surface changes by 
and thus the effect of the transverse shear is expressed by 


0 

a 
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(2.5) 


From (2.4) a compatibility equation is extracted in the form 






( 2 . 6 ) 


where e^^^ permutation symbol. 

When elastic shells are considered, Hooke's law yields 




(2.7) 


where h is the thickness of the shell. By use of the stress function 
F(X^) defined by 



p p F 

ay B6 ,y6 


( 2 . 8 ) 


the equilibrium equation (2.1) is satisfied, while (2.2) takes the 
form 


"^(xp.ap ^,ctP®ay®P<S*">Y'S ^ ^ ° 

The compatibility equation (2,6) becomes 


(2.9) 
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( 2 . 10 ) 


In general solutions of the system of differential equations (2.3), 
(2.9) and (2.10) will be quite complicated. However, considerable simp- 
lifications arise, if the differential operators can be factorized. Be- 
low only orthotropic materials with (X^) as principal axes are considered. 
This implies the following form of (2.7), 


"n 


r 1 

^1 

0 


Nil 

^12 

_ 1 
'F 

'’2 

‘^2 

1 

^2 

0 


^22 

2 Ci 2 


0 

0 

1 


^12 


( 2 . 11 ) 


where ” ^i^2’ 





In this case the factorization property amounts to 




( 2 . 12 ) 


Now introduce the geometrical mean values 
E = , V = /v-|V2 

and the orthotropy parameter 6 defined by 



In terms of these three parameters (2.11) takes the form 


(2.13) 


(2.14) 
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Following from an assumption of linear variation of the stress compon- 
ents a - over the thickness the moments are then given by 
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A linear relation between the angles and the shear forces V is 

0( Up 

assumed. In order to make elimination of the moments from (2.9) possi= 
ble this relation must be of the form 
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( 2 . 17 ) 


B is the effective transverse shear modulus, which in the isotropic case 

C 

may be taken as G [5]. 

The system of differential equations can now be simplified by elim- 
inating the moments. When the notation 
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Is adopted, the result is 


Vgf t nt 
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and 


(2.19) 
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( 2 . 20 ) 


with the extra conditions 
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3. vmEHSiomss parameters 

The following curvature measures are introduced 


1 _ 3=^Z 


R, 


3Xf ’ Rg 


1 _ 3^Z 
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3X| * Ri2 


SX^SXg 


(2.23) 


When dimensionless variables and parameters are defined as shown in 
Table 1, the differential equations {2.19)-(2.22) reduce to 


n4* 1 (yZ 3^ oiZ 3^ 




(3.2) 


v'*w + X^(l-KV^)(Xf - 2Xj^ + \l -^)$ = X‘*(1-kV^) ^ q 


(3.3) 
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and 

0-kV^)3x + K (-3^ - (3.4) 

0-KV^)3y + |^= K ^ (-3^ - -gjf ) (3.5) 

The usual static and geometric boundary conditions can also be 
formulated without explicit use of the orthotropy parameter 6. As a 
consequence the following presentation refers to an isotropic shell. 

The corresponding specially orthotropic solutions are easily found 
by use of Table 1. It should be pointed out that the shell parameter 
Xz is 6^ times the shell parameter used by Yuceoglu and Erdogan [12]. 

4. THE CVLIHVER - IWTEGRAL REPRESEmTiOW 

When considering the cylinder shown in Fig. 2, Xi=Xi2=0. The 
homogeneous equations corresponding to (3.2) and (3.3) then take the 
form 


v't - 0 = 0 

(4.1) 

v**w + (XX2 )Ni-kV2) 0 = 0 

(4.2) 

Elimination of either w or $ from (4.1) and (4.2) 
8th order differential equation 

leads to the same 

V®w + X^(I-kV^) = 0 

(4.3) 


Although the introduction of a finite transverse shear stiffness 
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represented by 1 /k does not increase the order of the differential 
equations (4.1) and (4.2), the order of (3.4) and (3.5) is increased 
thereby enabling the satisfaction of five boundary conditions as com- 
pared to four, when k=0. 

The procedure now is to give integral representations for ^(x.y), 
w(x,y), B^(x,y) and 3„(x,y). When a crack is present along part of 
the x-axis as shown in Fig. 2, different expressions must be given for 
the half-planes y>0 and y<0. Introduce the representations 



where i(5,y) and w(C,y) are given by 



(4.4) 

(4.5) 


(4.6) 


w(5.y) = j w(x,y)e”'^^dx 


(4.7) 


Due to (4.3) and the corresponding equation for $(x,y) both w(C,y) 
and 5(C,y) can be expressed as linear combinations of exp(mjy)* where 
m.=m.(?) are the roots of the characteristic equation 

J J 




(4.8) 


The notation 
p = m^ - 


(4.9) 
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Is now Introduced leading to the quartic equation 
P** + = 0 


(4.10) 


The solutions Pj(5), j=li..»4 to (4,10) are given 1n Appendix A. The 
Solutions nij(?) to (4.8) are selected such that 

Re[mj3<0 , “ -Wj , j=l,..,4 (4.11) 

When the displacement transform function w(C»y) Is given In the 


'* m.y 

I R,-(C)e ^ . y>0 

j=l ^ 

w(e.y) =<^ 8 m.y 

I R.UU ^ , y<o 

l3=5 ^ 


(4.12) 


the transform function 4(C*y) is found by substitution of (4.12) into 
(4.1). 

I ,x5.s ; ,f .2. , . “^y 


i(C.y) 


<t)' I . ^>0 

X jSl Pj 3 

■ <¥) I (f) Ri(«)« ^ • y«> 

j=5 Pj 3 


(4.13) 


Introduce the function H'(x,y) as 
38. 88 

w/v vl = ^ ^ 


The equations (3.4) and (3.5) give the equation 


(1- K ^ 72)'i' = 0 


from which the transform function ’f(?,y) is found to be 


(4.14) 


(4.15) 
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y(c,y) = 


Ai(5)e 
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[Agte)* 


. y>o 

‘‘2^ , y<0 


where 




2 + ^ 


k( i-vT 


(4.16) 


(4.17) 


Now ^(Xty) is assumed to be independent of $ and 0 . The equa- 

A y 

tions (3.4) and (3.5) are then solved, and the solution verified by 
substitution into the original equations. When a function x(>^*y) with 
the transform 

I ijf) Rj(5)e ^ , y>0 

x(5.y) y „ y (4.18) 

is introduced, the solution is found to be 


X 3x ^ 3y 


(4.19) 


3x 3x 


(4.20) 


The solution to the field equations has now been represented in 
terms of ten unknown functions A.(c), j=l,2 and R.(C)» j=l..-8, whici 

J J 

must be determined from the boundary conditions of the problem. 


5. BOumRv comrms - symmetric LOAvm 

For the axial crack shown in Fig. 2 the following five static 
quantities are prescribed at the crack surface. 
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The last expression is used in differentiated form in order to secure 
dimensional consistency. By substitution of the integral representa- 
tions of the ’‘•>' 5 *vi«)us section the following expressions are found for 
the half-plane y>0. The corresponding expressions for y<0 follow trivi- 
ally from changes of the indices. 

Nvv(x,y) = i I Pi'%(?)e J >e d? (5.2) 

1 r" X, 2 ** o m.y -i^x 

Vx.y) = ^ J (-iC=(f) E 

-00 J “ I 

-00 J“ I 

+ i?r^A.,(5)e’"’'^}e ^^^d? (5.4) 

M (x,y) =1^1 I ^ 

ioo J“ I 

+ ^[l+K(l-v)e']A.j(?)e Me d? (5.5) 
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V^(x,y) = ^ ^ 

1 -i?x 

+ •'^ ?^A-j(5)e /6 d^ (5.6) 

The problem of static boundary conditions on the crack In Fig. 2 
could In principle be formulated directly In terms of the unknown func- 
tions A.(c) and R.(C). This would lead to a system of coupled dual 

J J 

Integral equations, which is not easily solved even by numerical meth- 
ods. An alternative and numerically more suitable technique consists 
in the following. The Integral representations given In Section 4 are 
used to express geometric quantities, which are In a sense complemen- 
tary to the static quantities (5.T). The most direct choice Is the 
generalized displacements corresponding to the generalized forces (5.1). 

For the sake of brevity we shall limit our attention to the case 
of symmetric self-equilibrating loading. 

Nyy(x,y) = N^y(x>-y) * Nj^y(x»y) = -N^y(x,-y) 

Myy(x,y) = Myy(x>-y) » M^y(x.y) “ -Mjjy(x,-y) (5.7) 

Vy(x.y) = -Vy(x,-y) 

In this case 

Rj+4(?) = Rj(5) , A2(d = -A^(C) 


(5.8) 





{ 


whereby the number of unknown functions reduces to five. 


Due to (5.7) the boundary conditions for and are homo- 

Aj xy y 

geneous. We therefore only need to introduce generalized displacements 

* 

corresponding to N^y and Myy. By use of Hooke's law (2.15) and the 
equilibrium equations (2.1) we find 



i!i 

3y3 


= -lim 
y+O 


3X= 


(5.9) 


It is convenient to use the following two functions as unknowns, 


5**' “ yVo” 


(5.10) 


f(x) = lim 
y-»-0+ 


It ^y<’‘’J'> 


(5.11) 


Their integral representations are 



Rj(?)} e 

(5)+ 5"A^(0)e 



(5.12) 

(5.13) 


The five unknown functions A|(5) and Rj(?), j=l,..4 are now elim- 
inated by use of the three homogeneous boundary conditions and inver- 
sion of (5.12) and (5.13). It is noted that g(x)=f(x)=0 for x^;[-l,l]. 
After a few reductions we find 


A-,(?) 


-2 


fl it§ 
f(t)e dt 
il 


(5.14) 
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and the four equations 


I 5m.R.(5) = -1 ff(t)e^^^dt 

j=l J J 

u 

I P,.'^Cm-R.(5) = 0 
j=l J J J 

J. Pj ^fifnjRj(^) = i ^2^1 f(t)e^^^dt 


{5.15) 

(5.16) 

(5.17) 

(5.18) 


Due to the systematic nature of these equations the solution is 
straightforward. When use is made of the characteristic equation (4.10), 
we get 


*^1^^^ " em.(p^-P2)lp^-P3)(p|-:p4) ^“PifPi'^n-v)C^] 


^f(t)e'*^^dt 

-1 


+ (AXj^'^^l-icp,) |^g(t)e^‘^W (5.19) 

The expressions for R2(5), R3(0 and R^(?) are found by interchanging 
the indices. 

It should be noted that A-|(e;)= 0 does not follow from the limit 
process <+0. Special care must therefore be taken, when relating re- 
sults of the present theory to results from 8th order shallow shell 
theory. 


6, SINGULAR INTEGRAL EdUATJOMS 


T 





When A^(5) and Rj(?) from (5.14) and (5.19) are substituted Into 
(5.2) and (5.4), integral expressions for Nyy(x.y) and Myy(x.y) are 
obtained In terms oi the two unknown functions f(t) and g(t). For 
y=0 the integrals are defined by the limit y+0+. Contributions to the 
integrals, which are non-integrable for y=0 are extracted by use of 
the equations (5.14)-(5.18). 

The identity (5.17) is multiplied by 5[Cl (Az/X)^exp(-l5y[ ) and 
added to the integrand in (5.2). 


Vx,y) = ^ 

-CO J "“ ^ J 
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('"i+UDlyl -|5yl -1x5 

]}e e dC 


( 6 . 1 ) 


Asmymptotic expansion for large values of 5 yields 

,2 


|r| (ni4+Ul)ly| , 1 Pi 3 Pi X K‘*'l?niyl , , 

•)+i^e J ~ 1 - (1- --Oe J (6.2) 


Substitution of (6.2) into (6.1) and use of the identity (5.16) lead 
to a non-integrable term of the form 


Nyy ( X . O ) 


1 rco -Icyj i^(t-x) 

g(t) sgn(5)e e 
1 ^00 


-|5yl -1x5 

e d5 
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dCdt 


aM Ht 
t-x 


(6.3) 
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The last integral must be evaluated as the Cauchy principal value. 

By a similar, but slightly more complicated procedure, the follow- 
ing non-integrable contribution to the moment Myy(x,0) is found 


uS _ a 1-v^ 1 f(t) 

M^y{x,0) -^dt 


(6.4) 


After extraction of the non-integrable parts of the integrals the 
boundary conditions take the form of two singular integral equations. 

K^l(t-x)]g(t)dt + ^ 2TfNyy(x,0) 


-Kx<l 


(6.5) 


1 T 2 

I K 2 -,(t-x)g(t)dt + K 22 (t-x)]f(t)dt = 2ir ^ My^(x,0) 


-Kx<l 


( 6 . 6 ) 


When the symmetry properties m.(?)=m.(-?) and P^(5 )=Pj(-C) ai"e 

J J J J 

noted, the functions Kj|^{t-x) are easily brought in the following form 

- D^jlTT Prpjlr-aHPrP ?'^^ 

r p, + (i-v)e^ 

K^2(t-x) = K2^(t-x) = |q^2(-^) I — (p^-p2)(pY-'p3)(p|-P^) ^ 

sin?(t-x)d^ (6.8) 


K2g(t-x) = 


f“^X ^3 [p^ + (l-v)C^]^ 

1 (Pi-P2)(Pi"-^(PV^ 
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The summations imply interchange of the indices. 


Due to the presence of terms- of order the expression (6.9) is 
not suitable for numerical integration. The problem is solved by con- 
structing the following identity from (5.15), (5.18) and (5.19). 

+ 2k 11^ + 4 ^ = 0 

( 6 . 10 ) 

Kg 2 (t-x) can then be evaluated by the more suitable expression 


^ ' Y 'Pi(Pi“P2)(PrP3^^PrP4^ 


^22^ i'“X) 


f" K. t [pi + (l-v)?^] 

j (Pl-P2)(Pi“p3)"( Pi-P^ 


+ 2k ?(5+r^) + ■ ^ ■ ^''4^— ^ sinC(t-x)dC (6.11) 


In order to get continuity of the displacements outside the crack 
the solution must also satisfy the two conditions 


f g(t)dt = 0 
il 


A 

f(t)dt = 0 

"1 


( 6 . 12 ) 


The integral equations given above only apply when k> 0. In order 
to obtain the corresponding equations for the 8th order theory, where 
K=0, the function A-j(?) must explicitly be set equal to zero. Only 
M 's dependence on f(t) includes A-j(?), and the changes are therefore 
restricted to the last integral in (6.6). The coefficient (l-v^)A'* 
must be replaced by (l-v)(3+v)/X‘’, and the kernel K 22 (t-x) is now 


(6.13) 


KggCt-x 


) = j 12(f) X 




Cpi + (1-v)?^]^ ^ 

:ip^-P2)'(PTP3j(p^-^ 


+ } sine(t-x)d? 


7. THE ASy^PTOTIC STRESS TJELV 


The solution to the two singular integral equations (6.5) and (6.6) 
is of the form [13] 


g(t) = (i-t2)-'^G(t) 
f(t) = (l-t^)-^F(t) 


(7.1) 


In order to obtain the asymptotic stress field in a neighborhood of the 
crack tip x=^l the following formula, which is derived in Appendix B, is 
used 



{G(l)exp[i(C~sgn(5) J)] 


+ G(-l)exp[-i(^-sgn(?) J)]} 


(7.2) 


Substitution of (7.2) and the similar formula for f(t) into (5.14)- 
(5.18) makes possible an asymptotic analysis similar to the one which 
lead to the extraction of the Cauchy integrals. Now, however, (6.2) 
must be expressed in the slightly different form 


J j“"2 
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f 


1 


c. , j=l,2,.,, are constants, which do not influence the result. 

J 

The asymptotic expressions are obtained by application of the 
following formula [14; .3.944] 


("jM-l.-SS 

sin' 

io 

cos 
% / 


{6?)d? = 




sin 

cos 


(p arctg ~) 


(7.4) 


It is seen that a factor (?3)^ under the integral sign does not change 
the order of the result for 3,6-K). The order of the terms in (7.3) can 
then be evaluated as if |?y| were a constant. It turns out that only 
the first term of (7.2) contributes to the singular stress field around 
x=l . 


Substitution of (7.2) and (7.3) into (6.1) yields 


M (x v) ~ ^ 


^^(l+5|y|)e-5i^ls1n[5{l-x)-|]d5 


(7.6) 


In the same way we find 


N^^(x.y) ~ ^ ["i (i.5|y| ) e-5lj'lsin[S(l-x) - (7.6) 


0 


Nj^y(x,y) - 1^ ye“^l^lcos[^(l-x) - J]d5 


(7.7) 


The asymptotic expressions (7. 5)- (7. 7) are also valid for k= 0. For 
K>0 the asymptotic expressions for the moments are 


'’yy'’‘"^> ' ifelyz? 0+5|y|)e"'l^ls'"[£0-*) -fW? (7.8) 

Jo ^ ^ 
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M^j,(x.y) |"i {l.{|y|)e-«l^lsin[5(l-x)-J]d5 (7.9) 

V*’^> ■ llr 27^ lo ye"5l^'lcos[{(l-x) -f]d5 (7.10) 

It follows from i7.5)‘-{7 .^0) that the asymptotic membrane and 
bending stress fields are self similar for k> 0. This is in agreement 
with results reported in [1-4] and [11]. For k= 0 the asymptotic bend- 
ing stress field is found to be different and depend on Poisson's 
ratio V [7,10], 

When the coordinates r and 0 are defined by 

x-1 = r COS0 , y - r sin6 (7.11) 

application of (7.4) and use of trigonometric formulae lead to the 
result 

I cos(0/2) + \ cos{50/2) 

^cos(0/Z) - ^ cos(56/2) (7.12) 

I sin(0/2) + ^ sin(50/2) 

For specially orthotropic materials r and 0 are not polar coordinates. 
With the usual definition of the stress intensity factor 



we get 



*^ 22^^1 


(7.13) 
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= - ^ E [G(l) + zF(l)] 


(7.14) 


By inspection of the integral equations (6.5) and (6,6) it is seen that 
the stress intensity factor at the surfaces Z=+h/2 depends only on the 
parameters v, Xa and k. Thus the thickness effect and the influence of 
the effective transverse shear modulus are combined in the parameter 
k“E/(BX'’). 

8 , WUMERICAl SOLUnOW 


The singular integral equations (6.5) and (6.6) are solved by use 
of a quadrature formula of closed type developed in [15]. The integral 
equations are replaced by the algebraic equations 
n 




( 8 . 1 ) 


^22(‘r*k>]f(‘i» =f V'k.o* 

( 8 . 2 ) 


2h 


where 


t^ = cos (^ tt) , i=l,2,...,n 


^ 2 , w - 2 . * k— 1 j 2 j . • • »n~l 


and 


(8.3) 

(8.4) 


'^n*! '^n,n 2(n-l } 


_ 1 


n,i 


fpj- » i=2,3,. . . ,n-l 

-22- 


(8.5) 


i 


The extra conditions (6.12) are replaced by 

J, = “ (8-6) 

» 

The convergence of the method is estimated by evaluation of the 
coefficients bj in the following expansion 

G(t)-l^G(l)+^G(-l) + (l-t'‘) "f b.U.(t) (8.7) 

where Uj(t) is the Chebyshev polynomial of the second kind of degree j. 
The coefficients can be expressed explicitly in terms of the solu- 
tion to the algebraic equations [15], 

o n-1 1+t. 1-t- 

(e.8) 

The necessary values of the bounded kernels (t^-X|^) . 
and calculated by use of Filon's integration formula. 

9. RESULTS 

Numerical results are given in the form of stress intensity factors 
for two loading situations, constant membrane load and constant bending 
moment. For each loading situation two stress intensity factors are 
given, one for the average stress and one for the surface stress from 
bending. Introduce the following normalization of the stress intensity 
factors. For compressive membrane load 

= ( 9 . 1 ) 
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kb„ = [K,(h/2)-K,(0)](o^,^r)'’ (9.2) 

For a constant bending moment with maximum surface stress 

kbb = [Ki(h/2)-K,(0)](Ob /T)"’ (9.3) 

''rob " (9-’) 

The results are given as functions of the shell parameter X 2 

various values of a/h in Tables 2-5 and Figs. 3-6, The effective 

5 E 

transverse shear modulus •^f^has been used together with \)f=0.3. 

The first column of the tables and the dashed curve in the figures 
correspond to 8th order shell theory (k= 0). As expected from the dif- 
ference in the asymptotic moment fields the bending stresses show some 
differences. The membrane stresses from the 8th order theory, however, 
are found to be representative for very thin shells h/a<10. It is im- 
portant to note that in general 8th order shell theory gives non- 
conservative estimates of the membrane stresses. 

The extrapolated values of for are in good agreement 
with the results obtained in [2-4] for a plate. 
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APPENDIX A 


i 


\ 


Solution of the equation p** + (A25)''(1 “kp) = 0 
Following the method 1n [16] the solution of the quartic equation 
P** + (X25)’*0 -kp) = 0 (A.l) 

requires the real root of the cubic equation 

x" - 4{\2 ?)'*x - = 0 {A. 2) 

Now introduce the parameter 

n = 3/J i-f-) (A. 3) 

The real solution to (A. 2) is determined as the sum of 6i and S^ven 
by 

(A.4) 

This relation is rewritten in the fom 

3 e [exp (+ i Arccosn) . 0<ti<1 
®l,2 " ^ 

[exp (+ Arccoshn) » l^n 

The real root of {A. 2) is then given by 

fcos (j Arccosn) , 0<n<l 

X = 6, +02 = •] <A-6) 

[cosh (j Arccoshn) . l<n 
The roots p of (A.l) are found from 
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nI 


p^ + w^Tp + ix + ^ (X^S) 

^ C.VK 


1* _ 


= 0 


leading to 


P = ^.5 ft '(t)’ * /n7S o PT 6 73TJ 


(A. 7) 


0<ri<l n COS0 


(A. 8) 


For l<n cos( ) is replaced by cosh ( ). 
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APPENDIX B 


Asymptotic expansion of I(?) = j e^^^dt 


Rewrite the integral in the form 
1 

-1 
1 


I 


(5) = f ^ e’«*dt 


dt 


tU(t)+tj)(-t)]cos(5t)+iC(}»(t)-(|)(-t)]sin(5t)> 


Introduce the following series expansion of <j)(t) 

Mt) = I c.T.(t) 

0=0 ^ ^ 

and use the results [14; 7.355] 


|\2„^,(t)sin(5t) (-ir|02„,i(l5|)sgnU) 


I 


^T 2 „(t)cos(Ct) =(-1)2 J2nU5l) 


T.(t) is the Chebyshev polynomial of the first kind of degree j 

J 

J^(5) is the Bessel function of the first kind of order n. The 
then is 

1(5) = 'TT sgn(5)c2j+iJ2j+l^'^l^^ 

From the asymptotic formula [14] 


I 

( 


(B.l) 


(B.2) 


(B.3) 


and 

integral 

(B.4) 
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Vi^i> ” ASr<=o=(i5i-f--T)n+ou'’)] 


(B.5) 


we get 

1(5) ~/Hj^ J [CgjCosd^l -J) + i sgn(?)c2j.4.TSin(Ul -J)] (B.6) 

By use of T 2 p(+l)=l and "^20+1 leading term of the asymptotic 
expansion of 1(C) is found to be 

1(5) {<l>(l)exp[i(C-sgn{C) J)] 


+ <{i(-l)exp[-i(C“ sgn(C) J)]} 


(B.7) 


TABLE 1 


Dimensionless variables. 


1 

X - — ^ 


Coordinates 

X2 


^ * ? 


Displacements 


u - 

O 


6, ■ A 6, 


V . J-l!i 

Jl » 

®y '■i®2 


w 

^ a 


Stresses and stress resultants 


XX 


XX 


XX 


.^11 

JT 

0 

yy tT 

n .”’2 

°xy T 


6^22 
"yy ■ Tir 

N . "'2 
"xy fir 

- "11 
mz 

\y h*E 

M . ”'2 
"xy 


/5HB 


V^Vo 

“ “FT 

Parameters 




X- 

1 


2,** 






x;, - 12(.V) ^ 


12 


X- = 12(l-v*) ^ 
h 


1 E 
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a / h =2 


a / h =5 


, a / h =10 



0.0001 

0.0001 


0.0212 

0.0208 

0.0571 

0.0551 

0.0544 

0.0912 

0.0890 

0.0891 

0.1206 

0.1193 

0.1200 

0.1604 

0.1636 

0.1674 

0.1748 

0.1851 

0.1942 

0.1397 

0.1661 

0.1887 

0.0395 

0.0762 

0.1156 

- 0.1089 

- 0.0698 

- 0.0140 

- 0.2965 

- 0.2605 

- 0.1865 

- 0.4991 

- 0.4826 

- 0.3952 

- 

- 0,7369 

- 0.6343 



- 1.1829 


- 

- 









I 

i 

I 


Table 4. Stress intensity factors 


Xj 

K=0 

a/h=l 

— 

a/h=2 

a/h=5 

a/h=10 

0.01 

1.000 

0.747 


0.662 

mmm 

0.25 

0.998 

0.745 

0.698 

0.660 

mi 

0.50 

0.992 

0.738 

0.692 

0.655 


0.75 

0.983 

0.728 

0.684 

0.649 

mm 

1.00 

0.973 

0.716 

0.674 

0.641 

0.628 

1.5 

0.950 

0.693 

0.653 

0.623 

0.612 

2.0 

0.927 

0.671 

0.632 

0.605 

0.597 

3.0 

0.881 

0.633 

0.594 

0.572 

0.564 

4.0 

0.838 

0.603 

0-561 

0.540 

0.535 

5.0 

0.801 

0.580 

0.534 

0.513 

0.506 

6.0 

0.767 

0.562 

0.511 

0.489 

0.483 

7.0 

0.737 

- 

0.492 

0.468 

0.463 

8.0 

0.709 

- 

- 

0.450 

0.446 

10.0 

0.661 

- 

- 

- 

0.413 

12.0 

0.621 

- 

- 

- 

- 
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